Time Variation of the Fine Structure Constant Driven by Quintessence 
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There are indications from the study of quasar absorption spectra that the fine structure constant 
a may have been measurably smaller for redshifts z > 2. Analyses of other data ( 149 Sm fission rate 
for the Oklo natural reactor, variation of 187 Re /3-decay rate in meteorite studies, atomic clock 
measurements) which probe variations of a in the more recent past imply much smaller deviations 
from its present value. In this work we tie the variation of a to the evolution of the quintessence 
field proposed by Albrecht and Skordis, and show that agreement with all these data, as well as 
consistency with WMAP observations, can be achieved for a range of parameters. Some definite 
predictions follow for upcoming space missions searching for violations of the equivalence principle. 
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INTRODUCTION 

Independent observations of a number of absorption 
systems in the spectra of distant quasars (QSOs) seem 
to indicate that a, the fine structure constant of quan- 
tum electrodynamics, is slowly increasing over cosmolog- 
ical time scales P, Q . Specifically, the experiments indi- 
cate that averaged over redshifts 0.2 < z < 3.7, there is 
a 5.7cr deviation of the fine structure constant from its 
present value, namely Aa/a — — 0.57± 0.10 x 10~ 5 
On the other hand, terrestrial and solar system measure- 
ments provide several constraints on recent rate of vari- 
ation of a : (1) Analyses of the resonant fission reaction 
rate in the naturally occurring reactor at Oklo in Gabon 
provide a bound -0.9 x 10~ 7 < Aa/a < 1.2 x 10~ 7 
over the era z < 0.14 at the 95% CL 4j. (2) Based 
on plausible assumptions, new estimates of the age of 
iron meteorites (z ~ 0.45) combined with a measure- 
ment of the Os/Re ratio resulting from the radioactive 
decay 187 Re — > 187 0s have allowed a narrowing in the 
uncertainty of the average decay rate over the age of the 
meteorite This has been translated [f| into a strong 
bound Aa/a < 3 x 10~ 7 , following the original sugges- 
tion of Peebles and Dicke [jfl. (3) Recently, three years of 
observations of hyperfine spectra using atomic fountain 
clocks have allowed a la bound a/a < 1.6 x 10 -15 yr _1 
for the present time-rate of variation of a ||. (4) Ad- 
ditional bounds in this category have been derived ||, 
but these are weaker than the ones listed above. Fi- 
nally, there are constraints resulting from cosmological 
considerations: (1) Limits on the temperature fluctua- 
tions of the cosmic microwave background (CMB) could 
lead to a measurement with experimental sensitivity of 
\Aa/a\ < 10~ 2 - 10~ 3 , at z ~ 1000 Analysis of 

data from the Wilkinson Microwave Anisotropy Probe 
(WMAP)provides a bound -0.06 < Aa/a < 0.02 at 
95% CL \u$ i . (2) Big bang nucleosynthesis (BBN) con- 
siderations place bounds on |Aa/cc| on the same order of 
magnitude as those from CMB, though at much larger 
redshift, z ~ 10 9 - 10 10 



Over the last years, a second set of observations (most 
recently from WMAP has accumulated which indi- 
cate that the universe is spatially flat to within 1%. In 
addition, luminosity distance measurements of Type la 
supernovae strongly imply the presence of some unknown 
form of energy density, related to otherwise empty space, 
which appears to dominate the recent gravitational dy- 
namics of the universe and yields a stage of cosmic ac- 
celeration [l3 |. We still have no solid clues as to the 
nature of such dark energy (or perhaps more accurately 
dark pressure), but in recent years it has been associated 
with a dy namical scalar field cf> evolving in a potential 
V(<j>) [lif . generally called "quintessence" (lrij |. 

The universality of gravitational interactions implies 
that one may expect the Lagrangian below the Planck 
scale to contain non-renomalizable couplings of <f> to stan- 
dard model fields 0|- I n particular, the free Lagrangian 
for the electromagnetic field tensor will be modified 
to 



(1) 



with Mpi = (8ttG)~ 1/2 , the reduced Planck mass. On 
expansion about the present value 4>Q of <f>, this becomes 



r 

with A<p = 4> 
malization A, 
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<po and k = d^Zj, 
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The field renor- 



A^/Zp to obtain a canonical kinetic 

A/2 



energy, generates an effective charge e/Z F . Expansion 
to linear order about the present value eg, leads to 



Aa 
a 



_ A(f> 
'Mpi 



(3) 



Compatibility between the Oklo/metorite/atomic clock 
and QSO measurements can greatly constrain the dy- 
namics of <j>: its evolution should slow considerably be- 
tween the quasar era and the present epoch [lsj . 
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If the field 4> driving the variation in a is a quintessence 
field, then its evolution is further constrained by obser- 
vation. In particular, it must provide about 70% of the 
total energy density at present. Its equation of state, 
w<t> = P4>/P<j> (P4> = pressure, = energy density), 
is most strongly constrained by WMAP observations: 
W4, < -0.78 at the 95%CL 13j. Additionally, radiation 
dominance at the time of BBN must be maintained - no 
more than 20% of the energy density at that time can 
reside in quintessence 

In this work we search for a model of quintessence in 
which the evolution of the scalar field <f> in its potential 
V{4>) drives the variation in a [2£|. The particular case 
where <f> is the dilaton was examined in |2l|. With such 
a dynamics, however, the universal coupling of the dila- 
ton greatly constrains the variation in a, making it diffi- 
cult to comply with the QSO data. An alternative pro- 
posal [2^|. where the dynamics of the scalar field resides 
in a non-trivial Kahler potential, can provide a variation 
of a compatible with observation. However, in this work 
we wish to pursue a path based on sums of exponential 
potentials, which can be more explicitly linked to string 
theories (more on this below). It is also important to 
stress that our approach differs fundamentally from mod- 
els where the variations of the scalar field are prima rily 
driven by its couplings to non-relativistic matter l2jl 12 1 
(generically referred to as Bekenstein-type models 1230. 

Before proceeding, we take note of a discussion |26J 
which relates this type of variation in a to a large shift in 
the cosmological constant, and hence questions its viabil- 
ity. In the spirit of [23 we adopt here a more wait-and- 
see position, since present field-theoretic considerations 
all require fine-tuning to solve the cosmological constant 
problem, and may need to be totally supplanted (perhaps 
by self-tuning mechanisms, such as described in 28]). 

ALBRECHT-SKORDIS COSMOLOGY 
Quintessence Phenomenology 

An interesting model for quintessence has been pre- 
sented by Albrecht and Skordis |2^|. The scalar field 
evolves in a potential (hereafter we adopt natural units 
8ttG = 1) 

V{4>) = V p (<f>) exp[-A0] , (4) 

with an economic polynomial factor, 

V p (ct>) = (4>-B) +A, (5) 

in which the constants A,B,f3 and A are phenomeno- 
logically determined. Because of the polynomial factor, 
this potential differs in a critical manner from the much- 
studied pure exponential: although the tracking proper- 
ties are similar, it allows sufficient radiation dominance 



during BBN while evolving to quintessence dominance 
in the present epoch, all largely independent of initial 
conditions. For the particular case (3 = 2, the field is 
trapped in a minimum yielding a permanently acceler- 
ated universe if A 2 A < I, whereas if A 2 A > 1 (and for 
a small region 0<1 — A 2 ^4<Cl), the accelerated era is 
transient |30|. In what follows we set (3 = 2. 

Exponential potentials are ubiquitous in 4-D field the- 
ory descendants of string/M theory [3]]]. Additionally, 
sums of exponentials have been proposed |32l | in earlier 
attempts at stabilizing the dilaton and allowing super- 
symmetry breaking through gaugino condensates. In or- 
der to avoid the difficulties associated with identifying 
quintessence as the dilaton |2l| . in what follows we as- 
sociate 4> with moduli related to compactification, since, 
unlike dilatons, the latter need not be universally coupled 
to matter and gauge fields. The origin of the polyno- 
mial factor form of the potential in Eq. (0J can be linked 
to a non-trivial Kahler term in an effective 4-D super- 
gravity theory [33| . However, in such a case an explicit 
connection to string theory has not been successfully es- 
tablished. On the other hand, string/M theory does pro- 
vide motivation for superpotentials which are sums of 
exponentials [34J. A much-discussed [3^ recent exam- 
ple occurs in 1 1-dimensional supergravity, with a geom- 
etry consisting of a warped product of our 4 dimensional 
space-time and an internal compact 7-dimensional hyper- 
bolic manifold whose volume is proportional to the dual 
of the field strength. Upon dimensional reduction, the 
effective potential in 4 dimensions consists of two expo- 
nentials, the first proportional to the 7-dimensional Ricci 
scalar, the second to the volume of the compact space. 
The exponents are proportional to the breathing modu- 
lus of the warp factor. 

With this in mind, we can think of the Albrecht- 
Skordis potential as the limiting case of three exponen- 
tials 



V(4>) = Ce' (A+e)( ^ B » -2C(I-e 2 ^/2)e- A ^- B ) 

+ Ce -(A-e)(0-S) j (6) 



where e 2 C = e~ AB , and e 2 <C 1. In order that the re- 
duction to a 4-D theory generate three exponentials with 
a single modulus, it is necessary to split the compact- 
ification of the internal space so that the resulting in- 
ternal manifolds have different warping factors with a 
single collective coordinate. Additional dilatonic degrees 
of freedom can arise if one starts from I0-D string the- 
ory (such as Type HA or Type IIB), and these need to 
be constrained or stabilized. We are currently studying 
these possibilities and the results will be presented else- 
where 361. 
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Cosmological Evolution of a 



and 



In order to follow the evolution of the fine structure 
constant, we need to trace the temporal behavior of <p 
since the quasar epoch, while at the same time requiring 
that the field provide a successful model for quintessence. 
This evolution follows from the 4-D effective action 



5' 



a x \l—q { — 



-g\- + C <j> + C vis + [Z F (<j>) - 1] C c 



(7) 

As usual, R is the Ricci scalar derived from the metric 
tensor g^, 



1 



(8) 



£ V i S is the Lagrangian of visible matter (baryons, pho- 
tons, and also baryonic and neutrino dark matter), and 



I p p^ v 



(9) 



Note that we have omitted terms coupling (j> to the 
other matter fields because we are searching for a non- 
Bekenstein solution, i.e., the potential is primarily driv- 
ing the evolution of </>. 

The equation of motion for then reads 



3Hd> = 



dV 



where 



if = I - 

a 



- <j) 2 + V(<f>) + pm + Prad 



(10) 



(11) 



with p m and p r the matter and radiation energy densi- 
ties, respectively. It is more convenient to consider the 
evolution in 



x = In a = — ln(l + z) , 



(12) 



with the present value of the scale parameter oq = 1. 
Denoting by a prime derivatives with respect to x, the 
equation of motion for (f> becomes 



pi # /2 +3dtV 

+ 3</>H — = 0, (13) 

V + p 



1 - ^76 



where p = p m + p T . Quantities of importance are the dark 
energy density 



1 



P4> 



H 



2 J/2 



+ v, 



(14) 



generally expressed in units of the critical density (fl 

P/Pc) 



3H 2 



(15) 



V + p 



The equation of state is 
' H 2 d> 12 



V 



H 



2 JJ2 



V 



(16) 



(17) 



For p m and p T we adopt the expressions 

p m =Ce- 3x , p r = 10- 4 Ce- 4a: /(x), (18) 

where C = Q m ,oPc,o, and f{x) parameterizes the x- 
dependent number of radiation degrees of freedom. In or- 
der to interpolate the various thresholds appearing prior 
to recombination (among others, QCD and electroweak) , 
we adopt a convenient phenomenological form, 



f(x) = exp(-a;/15). 



(19) 



In natural units p c .o — 2.3 x 10 120 Hq. Finally, we set 
our constants O mi o = 0.3, h$ — 0.7 in accordance with 
WMAP observatio ns [H . 

As described in |29J, for a wide range of parameters 
this potential allows a plausible cosmological behavior 
independent of initial conditions. For definiteness, in this 
work we take the initial kinetic energy in the field <f> to 
equal its initial potential energy. The remaining degrees 
of freedom A and A will be used in order to study the 
variation of the fine structure constant over the history 
of the universe. We have found that the recent variation 
of 4> is most directly controlled by the curvature of the 
potential at its minimum. This is given by the mass of 
the scalar field, 



2 (1-/C)e 



-XB—K 



where 



K = l- VI - AX 2 



We will state our results in terms of 



-2 _ 



m; 3n w A 2 (1-/C) 



Hi 



K 



(20) 



(21) 



(22) 



For fixed values of m 2 and A, the value of B is fixed by 
requiring that the current dark energy density constitute 
a fraction fi^o of the critical density. 

Now, Eq. I|13|l is integrated for a range of values of A 



and 



frc 



- io- 30 



to the present epoch. Of these, a 



small subset has been found to be of interest with respect 
to the data on the variation of a presented in the Intro- 
duction. Within this subset, two types of solution can 
be identified: in the first, the variation of a shows an os- 
cillatory behavior between the present and QSO epochs, 
which reflects (through Eq. (J3J) the oscillatory behavior 




FIG. 1: The solid line indicates the variation of Aa/a for 
m% = 52.5, A = 8.5, B = 32.0, and « = -3.3 x 10 -4 . The 
cosmological evolution of a is superimposed over a binned- 
data sample from 72 QSO absorption systems: The points 
indicated by * correspond to two HI 21 cm and molecular 
absorption systems |37j |. Those points assume no change on 
the proton p-factor, and therefore should be interpreted with 
caution. The 7 squares are binned results for 49 QSO absorp- 
tion systems |38|. The lower redshift points (below z ~ 1.6) 
are based on (Mgll/Fell) and the higher redhsift points on 
(Znll, CrII, Nill, A1III, A1II, Sill). The single point indicated 
by • represents the average over 21 QSO SilV absorption 
doublets using the alkali doublet method |39|. The embed- 
ded box details the behavior of Aa/a for small values of z. 
The solid horizontal lines indicate the bound derived from 
the nuclear reactor at Oklo, whereas the dotted line in the 
middle indicates the predicted average value of Aa/a within 
< z < 0.14. 



of in the potential well; in the second, the motion of 
4> as it comes to its present value is overdamped, so that 
the change in a is a monotonically decreasing function of 
z. 

In Fig. ^ we show an example of the oscillatory be- 
havior, for the set of parameters m| = 52.5, A = 8.5, 
B = 32.0, and k = -3.3 x 1CT 4 . Although the fit to the 
QSO data is not particularly good (x 2 /d.o.f. = 16/7), 
a solution of this type (with some V p of higher order in 
Eq. can lead to an interesting prediction for future 
observations: |Aa/a| diminishes with increasing redshift 
beyond z = 3. The fit is consistent with the Oklo bound 
(see the inset), but it fails to comply with the new mete- 
orite bound. 

Further comparison with data must include compli- 
ance with measured bounds on violation of the equiva- 



FIG. 2: The upper panel shows the evolution of <p for m| = 
52.5, A = 8.5, B = 32.0 (solid line) and fh% = 10.5, A = 
6, B = 45.4 (dotted line) as a function of log 10 (a). Today 
corresponds to a — 1, for BBN a w 10~ 10 , and for the Planck 
epoch a « 10 -30 . The second panel shows the evolution of £1$ 
for the ^-parameters described above. The lower panel shows 
the evolution of the equation of state wj, (same conventions 
than the upper panel). These two solutions are consistent 
with the BBN requirement Q.^, < 0.2 [ToT ]. show the established 
radiation and matter dominated epochs, and at the end yield 
an accelerated quintessence era. 



lence principle, in the form of limits on composition de- 
pendent inertial forces [2lL I23I l4(ij . These can translate 
into an u ppe r bound on k. Along these lines, Olive and 
Pospelov |23J parametrize the ^-dependent modification 
of the electromagnetic Lagrangian via an equation anal- 
ogous to Eq. J2J, 

£ cm = -\ (1 + Cf^ + • • •) F^F^ , (23) 

where M* is the analogue of Mp\ in the </> sector, and 
the field <j> in Eq. (|23|l is defined to have a canonical 
kinetic energy. Comparing Eqs. @ and (|23|l . we find 
k =_Cf/v / 2w, where cu = M%/2M£ l . The limit derived 
in |23J Cf/V^ < 10~ 3 , requires 

k<7 x 10~ 4 . (24) 

We can see that our value, k = 3.3 x 10~ 4 , is con- 
sistent with the one given in Eq. (|24J) . but close: tak- 
ing this model seriously would suggest that composition- 
dependent inertial forces could be observed in an im- 
proved round of experiments. 
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The third constraint mentioned in the Introduction is 
the variation of atomic clock measurements. In terms of 
our variables, the fractional change in the present epoch 
is 



— = kH 4> , 

a 



(25) 



where for this set of parameters 



-16 



yr 



>' = -7.8 x 1(T 3 . This 
well within the bound 



implies a/a = 1.8 x 10 
stated in the Introduction. 

Finally, as can be seen in Fig. the change in </> be- 
tween BBN and the present is ~ 10. On the assumption 
that the dominant variation in Zp((f>) is linear over this 
domain, the fractional change in a is well within the sen- 
sitivity of present analysis. A similar statement holds for 
the CMB lU. 

We now turn to our second example, shown in Fig. [3] 
in which rh\ = 10.5, A = 6, B = 45.4, and k = 2.1 X 10 -4 . 
This solution passes muster on several counts: (1) The 
fit to the QSO data is acceptable - x 2 /d.o.f — 9.8/7, 
corresponding to 20% CL. (2) As can be seen in the em- 
bedded box in Fig.0 the low-z model results are in good 
agreement with both Oklo and meteorite constraints. It 
is also consistent with equivalence principle bounds, the 
value of k being more than a factor of three below the lim- 
its discussed above. It should be stressed that less than 
an order of magnitude improvement in the experimental 
sensitivity for composition- dependent inertial forces may 
provide a direct test for this type of model. As in the 
previous case, this solution is also consistent with limits 
imposed by atomic clocks: we find 4>' = —2.4 x 10~ 3 , 
giving a/a — 3.5 x 10~ 17 yr -1 . (3) Finally, the model 
also gives variations in a during the BBN and recombi- 
nation epochs which are well within present sensitivities: 
for z = 1100 (CMB), we find (retaining the linear expan- 
sion of Z F (4>)) that Aa/a = -4.4 x 10~ 4 ; for z = 10 10 
(BBN), we obtain Aa/a = -0.0025. 



DISCUSSION AND CONCLUSIONS 

(1) In this work we have attempted to associate the pos- 
sible temporal variation of a, the electromagnetic fine 
structure constant, as indicated in absorption spectra of 
QSOs, with the temporal evolution of the quintessence 
field which is responsible for the present dark energy con- 
tent of the visible universe. The outstanding obstacle in 
maintaining this association is the observation that Aa, 
the deviation of a from its present value, is much smaller 
in the recent past (say z < 0.45) compared to the varia- 
tion indicated by the higher- z QSO data. The reconcili- 
ation can be effected if the quintessence field has under- 
gone a rapid slowing in the recent past. The quintessence 
model that we study (the Albrecht-Skordis model) has 
precisely such a property. We illustrate our results with 
an example (Fig.[3J) which can simultaneously fit the QSO 
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FIG. 3: Cosmological evolution of Aa/a for = 10.5, A = 
6, B = 45.4, and k = 2.1 x 10~ 4 . The experimental data 
points are those of Fig. The embedded box displays both 
bounds from the Oklo reactor and from meteorite analyses 
(horizontal solid lines). The predicted average value of Aa/a 
within z < 0.14 and z < 0.45 is indicated with dotted lines. 



data and comply with the upper bounds on Aa from the 
Oklo and meteorite analyses and atomic clock measure- 
ments. This class of solutions presents several inevitable 
predictions: (a) deviations from universal free fall should 
be observed when experimental sensitivity is improved 
by a factor or 10 (6) QSO measurements at ever-larger 
redshifts should continue to show a monotonic decrease 
in Aa : in this type of solution the quintessence field is 
just receding from its first turning point in its damped 
oscillation about the fixed point. Another type of behav- 
ior (example in Fig. ^) is seen in solutions which do not 
satisfy the meteorite bound: there the quintessence field 
has completed a few oscillations, allowing a return of a 
to its present value at large redshift. 

(2) The evolution of the quintessence field, and there- 
fore of a, is determined by the quintessence potential. As 
noted in the text, this is in contrast to Bekenstein-type 
models where changes in a are driven by the matter den- 
sity. The latter allows local spatial variation |42][. which 
could provide different dynamics for Aa in our local en- 
vironment (Oklo and meteorites) and over cosmological 
scales (QSO data). In the quintessence model, only the 
zero mode is relevant, and this option of spatial variation 
is much suppressed. 

(3) Comparison with Ref. [24^ . in which the evolution 
of 4> is driven by both matter and by a harmonic scalar 
potential, shows that we require a much larger value of 
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fn^ for agreement with data. The origin of this difference 
is of interest, since it highlights the constraints imposed 
by requiring that <f> fulfill its role as a quintessence field. 
Small values of the mass (< H) will flatten the potential 
well to such an extent that the field escapes entrapment 
and leads to an exit from the de Sitter phase. Moreover, 
in our model, cj> is not initially placed at the equilibrium 
point, and thus for a shallow well will have considerable 
velocity during the present era, leading to a strong dis- 
agreement with the low z data. 

(4) Because the coupling k of the scalar field to the 
electromagnetic Langrangian plays a dual role in deter- 
mining both the variation of a and the violation of the 
equivalence principle, planned experiments on univer- 
sal free fall can directly test the viability of our model. 
These include the mission MICROSCOPE from the Cen- 
tre National d'Etudes Spatiales (CNES) expected to fly 
in 2005 0, and the National Aeronautics and Space 
Agency (NASA) and European Space Agency (ESA) 
mission STEP (Satellite Test of the Equivalence Prin- 
ciple) 0. 

(5) Our prediction for the present rate of variation of 
a is an order of magnitude below present atomic clock 
sensitivities. Improvements in the accuracy of such mea- 
surements to the range of 10~ 16 yr _1 , of the order of our 
result, are anticipated for the near future ||. 

(6) Amusingly, we can project the evolution of a forward 
in time. For our model in Fig. [31 <j> is at a turning point in 
the potential well of V and consequently <f> (and a) would 
reverse motion in the future. For the model in Fig. cf> 
is already decreasing at present and will continue to do 
so in the near future. 
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